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Abstract 

An empirical model has been developed for predict- 
ing the temperature distribution downstream of a row of 
cool jets injected normal to a hot confined crossflow 
The model is based on the assumption that all properly 
non-dimensionalized vertical temperature profiles can 
be expressed in a self-similar form. The scaling 
parameters in this form have been correlated in terms 
of the independent flow and geometric variables. The 
effect of parametric variation of each of the independent 
variables on the experimental and predicted profiles are 
examined. The predicted distributions show excellent 
agreement with the data over a wide range of the inde- 
pendent variables. 

Introduction 

The objective of this study was to model the pene- 
tration and mixing characteristics of multiple jets of 
cooling air injected normally into a heated crossflow in a 
constant area duct. The study was motivated by consid- 
erations of dilution zone mixing in gas turbine combus- 
tion chambers. For this application, rapid mixing of the 
diluent air with the hot combustion gases leaving the pri- 
mary zone is desired to (1) provide a rapid quench for 
any continuing chemical reactions. (2) provide a suitable 
temperature distribution at the turbine inlet, and (3) re- 
duce combustor length. 

The experimental data on which the model is based 
were obtained by the Aerojet Liquid Rocket Company 
(ALRC) under NASA Contract NAS3-15703, and are re- 
ported in Refs. 1 and ?. The temperature field model 
was developed by ALRC under NASA Contract NAS3- 
18026 (Ref. 3). Data for selected tests from the Ref. 2 
experiments have been used by Pratt & Whitney Aircraft 
to generate correlations to characterize the behavior of 
the dilution process (Ref. 4). The results of Refs. 3 
and 4 are similar in that they both provide satisfactory 
profi-u predictions over the range of the data base used 
in Ref. 4. Because a larger data base was used in Ref. 
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3, these correlations can be used for a wider range of 
independent variables. Also, the correlations of Ref. 3 
are somewhat less cumbersome than those in Ref. 4, 
and the relative importance of each of the independent 
variables in the correlations is more readily appaient. 

In this paper, the effect of parametric variation of each 
of the independent variables on both the experiments 1 
and predicted profiles (using *he correlations of Ref 3) 
are examined. 

The Data Base 

A schematic of the multiple jet flow field with the 
principal flow and geometric variables identified appears 
in Fig. 1. The test section was 30.48 cm wide by 
10. 16 cm high. For all of the test results shown in this 
paper the nominal mainstream conditions were: velocity, 
U — = 15 m/sec; temperature, T^ * 600 K. The Jet 
velocity, Vj, was varied from 25 to 85 m/sec with the 
jets at ambient temperature, T. = 300 K. The Jets en- 
tered the test section through snarp-edged orifices in the 
plate separating the main duct from the secondary air 
plenum chamber. 

The test conditions were established with the jet-to- 
mainstream momentum flux ratio, J » (pjV^/jp^U^), 
and the jet-to-main stream density ratio, Pj/p^ , as the 
primary independent flow variables. Momentum flux 
ratios investigated ranged from 6 to 60, with density 
ratios varied from 1. 5 to 2. 5. The symbols and param- 
eters referenced in this paper are defined in Appendix A. 

The primary independent geometric variables were 
the orifice size and spacing between adjacent orifices. 
These were expressed in dimensionless form as the ratio 
of the duct height to the orifice diameter, (4 < H/D < 16), 
and the ratio of the orifice spacing to the orifice diam- 
eter, (2 < S/D < 6). A though the experiments were per- 
formed with S/D and H/D as the independent geometric 
variables, the ratio of these, S/H, may be considered as 
an independent variable in the ^ace of either S/D or 
H/D. The orifice area may be expressed in dimension- 
less form as Aj/A^ * (ir/4)/[ (H/D> (S/D^ J which is the 
orlfice-to-mainstream area ratio. The orifice configu- 
rations investigated are shown in Fig. 2. The ratio of 
the Jet flow to mainstream flow, Wj/w^ , can be ex- 
pressed in terms of the independent flow and geometric 




variables as 



For this investigation, the orifice discharge coefficient, 
C d , varied from n.66 at the lowest momentum flux ratio 
to 0.62 at the highest momentum flux ratio. 

Because the objective of the study was to identify 
orifice configurations appropriate for optimum mixing 
within a minimum combustor length, the downstrerm 
stations surveyed were defined in terms of the duct 
height, H. Measurements of total pressure and temper- 
ature were made at 20 vertical and 21 horizontal posi- 
tions in each of five planes at downstream locations from 
X/H * 0 25 to X/H - 2.0. The total span of the 21 hor- 
izontal positions was varied to correspond to twice the 
orifice spacing for each plate. This distance varied 
from 2. 54 cm for the smallest orifice spacing to 20 32 
cm for the largest spacing. The Z plane through the 
orifice center is defined as the centerplane (Z ■ 0) and 
the Z plane midway between adjacent orifices is defined 
as the midplane (Z - S/2) . Thus the span of the meas- 
urements provided 4-fold data redundancy. Additional 
aspects of the experimental program and the facility 
used are discussed in Ref. 2. 

The results for the temperature field are presented 
as vertical profiles of the dimensionless temperature 
difference ratio, 0, where 


In this expression, 0 is the local temperature difference 
ratio given by Eq. (2), and 0 c , 0^^, Y^H, and 
*1/2 /H are scaling parameters as shown in Fig. 3. 0 p 
is the maximum temperature difference ratio in the ver- 
tical profile, and Y c is its location. The line defined 
by the locus of Y c as a function of downstream distance 
(X) for Z»0 is the thermal trajectory (centerline). 

Because the flow is confined and the profiles are not 
symmetric about the centerline, the half-widths (W* , 9 ) 

and the minimum temperature difference ratios (0^^) 
are different for the + side (Y/H > Y^H) and the 
- side (Y/H < Y^H) of the profiles. Correlations have 
been developed for each of the scaling parameters in 
terms of the independent variables J, S/D, H/D, X/H, 
and Z/S. The complete set of correlations are given in 
Appendix B. This is only one possible set of correla- 
tions; in Ref. 4, Cox gives another set which also yield 
satisfactory profile predictions. The primary differen- 
ces between the results given here and those in Ref. 4 
are (1) the forms of the present correlations are simpler 
than those of Ref. 4 and the relative importance of each 
of the independent variables is more readily apparent 
and (2) profiles can be predicted using the correlations in 
Appendix B over a wider range of the independent varia- 
bles since the data base used in the present investigation 
was larger than the one used in Ref. 4. The main weak- 
ness of the present correlations is that their form pre- 
cludes their use in predicting profiles for semi-confined 
flows (large H/D or large S/D) or the single jet flow 
(large H/D and S/D). 
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Results and Discussion 


Effect of Downstream and Lateral Distance 


and T is the local total temperature. T^ is the total 
temperature of the undisturbed mainstream flow', and Tj 
is the Jet total temperature. Because T^ Tj, the 
largest values of 0 in any profile correspond to the 
coolest regions of the flow. 

Flow Field Model 


The empirical model for this three-dimensional 
shear flow is based on the observation that the properly 
non-dlmensionallzed vertical temperature profiles every- 
where in the flow' field can be expressed in the following 
self-similar form: 



The variation of the centerplane (Z = 0) temperature 
profiles with downstream distance tor both the smallest 
(J * 6) and largest (J = 60) momentum flux ratios tested 
are shown in Fig 4 for an orifice configuration with 
S/D “ 4 and H/D =8. For this geometry, the agreement 
between predicted and experimental profiles is slightly 
better at the low momentum flux ratio than at the high 
one. At the high momentum flux ratio, the jet penetra- 
tion increases ir.onotonically with increasing downstream 
distance. However, at the low momentum flux ratio, the 
penetration is non-monotonic, with a slight recurving 
toward the injection wall evident at the largest down- 
stream distance. This effect appeared whenever the jet 
penetrated to less than half the duct height, and vas par- 
ticularly acute when a two-dimensional jet flow was es- 
tablished near the injection wall (Y * 0). In Fig. 5. la- 
teral profile 8 at X/H * 1 are shown for the same condi- 
tions as in Fig. 4. For the low' momentum flux ratio, 
the profiles vary substantially with lateral distance, and 
the moJel provides a satisfactory prediction of both the 
decrease in the maximum penetration and the decrease 


in the maximum temperature difference with increasing 
distance from the centerplane. A thermal contour plot 
for this condition would show the traditional "kidney" 
shape resulting from the twin vortices, which develop 
with the flow from each Jet. The lateral variation of the 
temperature profiles is small at the high momentum flux 
ratio, indicating that there is substantial interaction be- 
tween adjacent Jets, which tends to diffuse the vorticity. 
The result is a ne: ' ly plane flow as shown in Fig. 5(b) . 

Effect of Flow Variables 

The jet-to-mainstream momentum flux ratio, 

J -<P,vf)/<p«l£>. was found to be the most important 
operating variable influencing the mixing. The magni- 
tude of this effect on the centerplane (Z - 0) profiles is 
shown in Fig. 6. These profiles are at a distance of one 
duct height downstream from the injection point for a 
geometry with S/D = 4 and H/D * 12 . The Prediction 
of the centerplane profiles is good except that the + side 
half-width is too large for low J values. It should be 
noted that the ratio of jet flow to mainstream flow in- 
creases by slightly more than a factor of three from 
J = 6 to J = 60. 

Although experiments were performed with the den- 
sity ratio varied from 1. 5 to 2. 5 at constant momentum 
flux ratio, no consistent correlation with jensity ratio 
was found in the present investigation. Thus the density 
ratio does not appear in any of the co relations in Appen- 
dL; B. This is not meant to imply tb it the mixing is in- 
dependent of density ratio for all Pj/p* , but rather in- 
dicates that, for the limited variation in the experiments 
(1. 5 < (pj/p^) <2.5), the effect, if any, is of second 
order importance. 

Effect of Geometric Variables 

The momentum flux ratio, and the density ratio are 
dependent on the combustor design conditions, e g. , 
inlet temperature and pressure, reference velocity, 
pressure drop, and fuel-air ratio. Thus for a given 
combustor, the dilution zone design parameters which 
may be varied to influence the combustor exit tempera- 
ture distribution are the orifice diameter, orifice spac- 
ing, and orifice shape. However, since the flow condi- 
tions for a given combustor determine the required ori- 
fice area, the orifice size and spacing must be correctly 
coupled. In Fig. 7, centerplane temperature profile 
data at X/H - 1 for momentum flux ratios of 6 and 60 
are shown for three orifice configurations all having 
oriflee-to-mainstream area ratios of 0.049. The differ- 
ences in the profiles are evident, indicating that for a 
given operating condition considerable variations in the 
downstream temperature distributions can be effected by 
changes in the orifice diameter and spacing. The most 
significant difference between the experimental and pre- 


dicted profiles in Fig. 7 is for the largest orifices at the 
highest momentum flux ratio. These conditions cause 
flow impingement on the opposite wall, and the tempera- 
ture distributions in the resultant strongly bifurcated 
flow do not conform to the self-similar profile shape as- 
sumed in the model. Thus, the lack of agreement be- 
tween predicted and experimental profiles is not sur- 
prising. 

The variations shown in Fig. 7 are consistent with 
the results of Ref. 5, where fewer large orifices were 
found to provide better Jet penetration than a larger num- 
ber of small orifices for the same orifice area. In the 
following paragraphs, the effects of parametric variation 
of the geometric variables at a constant momentum flux 
ratio of J * 25 are discussed. 

Effect of varying orifice spacing at constant orifice 
diameter (H/D = 8. J = 25, X/H - 1) - For a given Jet 
diameter, momentum flux ratio, and downstream dis- 
tance the penetration of the dilution Jets increases with 
increasing spacing as shown in the centerplane profiles 
in Fig. 8(a). Also, the maximum temperature differ- 
ence in each profile decreases as would be expected 
since the ratio of Jet flow to mainstream flow varies in- 
versely with S/D. The agreement between experimental 
and predicted profiles is excellent in all cases. The 
variation of the corresponding midplane profiles is shown 
in Fig. 8(b). As spacing increases, the centerplane and 
midplane profiles differ as the flow develops three- 
dimensionally . The prediction of the midplane profiles 
is good, but of courc? there is some spacing for which 
further increases woulc have no effect on the experi- 
mental distributions. Beyond this spacing, the flow is 
semi-confined and the model (as presently formulated) is 
not appropriate. 

Effect of varying orifice diameter at constant 
S/D(S/D =2, J = 25) - For a given S/D, momentum flux 
ratio, and downstream distance, the Jet penetration in- 
creases significantly with increasing orifice diameter 
(decreasing H/D) as shown in the centerplane profiles at 
X/H * 1 in F*g. 9 The maximum temperature differ- 
ence in these profiles increases slightly as the ratio of 
Jet flow to mainstream flow increases (Wj/w^ ~ l/(H/D)), 
but the primary effect is in the percentage of the span 
cooled by the dilution flow. The predicted and experi- 
mental profiles compare favorably except for the small- 
est orifice size where the model underpredicts the tra- 
jectory (too much recurving) and underpredicts the maxi- 
mum temperature difference. 

Some insight into the confining effect of the opposite 
wall is given by the centerplane profiles of 0 vs Y/D 
at constant X/D (Fig. 10) The maximum temperature 
difference decreases slightly with decreasing H/D, but 
the primary effect is the suppression of the Jet penetra- 
tion. The effects shown here are in agreement with the 
results of experiments on heated Jets in a confined cross- 
flow reported in Ref. €. 
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Effect of varying orifice diameter at constant 
spacing (S/H = 0.5. J = r :5) - For a constant S/H (con- 
stant spacing), constant momentum flux ratio, and con- 
stant downstream distance, increasing the orifice diam- 
eter shifts the temperature profiles to higher 0 value? 
consistent with the greater dilution airflow. Although 
the penetration of the jets also increases slightly, the 
shape of the temperature profiles is not altered appre- 
ciably. This effect is shown by the predicted and <jcperi- 
mental center-plane profiles at X/H = 1 and 2 in Fig. 11 
for three orifice configurations with a constant 
S/H = 0. 5. but with jet flow- to-mainst ream flow ratios 
varying from 0.11 to 0.46. The profiles shown in Figs. 

8, 9, and 11 suggest that for J - 25, optimum mixing of 
the jets with the mainstream flow is achieved with 
S/H = 0.5. The profiles in Fig. 7 show that for J * 6, 
the best mixing was achieved with S/H * 1, and fer 
J = 60, the best mixing was achieved with S/H = 0.25. 

Application to Combustor Design 

. These results suggest that for a given momentum 
flux ratio and downstream distance, combustor design 
procedure should first identify the S/H value required 
to obtain the desired profile shape, and to position the 
profile at the desired spanwisc location. The orifice 
size would then ue chosen *o provide the required jet- to- 
mainst ream flow ratio. Because the penetration varies 
slightly with orifice size, adjustments (using the corre- 
lation relations in the flow model) would be necessary to 
arrive at the final dilution jet configuration. Although 
the geometries and test conditions of the experimental 
data base are representative of current design practice 
for annular gas turbine combustors, the flow model can- 
not completely describe dilution zone performance, 
since, as discussed by Cox in Ref. 4, the effects of liner 
cooling airflow, non-uniform dilution zone inlet tempera- 
ture distribution, and flow' area convergence in the dilu- 
tion zone are not considered. Also, die experiments 
considered injection from one side toward an opposite 
wall, thuu the results are applicable to one-side entry 
''ombustors or (since the wall may be considered to be a 
plane of symmetry as shown by the opposed Jet and Jet/ 
wall data in Ref. 6) to dual-side entry combustors wiJi 
directly opposed Jets. At present there arc no experi- 
mental data available that could be used to extend the 
model to dual-side, staggered jet dilution zone configu- 
rations. 


Summary of Results 


An empirical model is presented for predicting the 
temperature distribution downstream of a row' of dilution 
lets injected normal to a hot confined crossflow. The 
model is based on the assumption that all properly non- 
dimen sionalized vertical temperature profiles can be 
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expressed in a self-similar form. The scaling param- 
eters in this form have been correlated in terms of the 
independent variables. These are the momentum flux 
ratio (6 < J < 60) , the ratio of the spacing betwe® adja- 
cent orifice to the orifice diameter (2 < S/D < 6) , the 
ratio of the duct height to the orifice diameter 
(4 < H/D < 16), the dimensionless downstream distance 
(0.25 X/H < 2), and the dime 'S ionless off-ccnterplane 
distance (0 < Z/S ^ 0. 5). Profiles predicted using the 
flow model are in excellent agreement wrlth the experi- 
mental data except for combinations of the flow and geo- 
metric variables which result in strong Jet impingement 
on the opposite wall. 

The effects of parametric variation of the indepen- 
dent variables on the predicted and experimental profiles 
can be summarized as follows: 

1. The momentum flux ratio was the most important 
flow variable influencing the penetrations and mixing. 

For the limited range of density ratios covered by the 
experiments, no correlation with density ratio was found. 

2. For a given orifice diameter, increasing the 
spacing between adjacent jets increased the penetration 
and increased the uniformity of the vertical profiles. 
However, horizontal distributions oecame more non- 
uniform with increased spacing as the mixing became 
more three-dimensional. 

3. For a constant spacing to diameter ratio (S/D), 
penetration and mixing at any X/H increased with in- 
creasing jet diameter. 

4. For a constant spacing (S/H), increasing the ori- 
fice diameter increased the magnitude of the tempera- 
ture differoice ratios, and increased the penetration 
slightly, but the basic shape of the profiles was not al- 
tered appreciably. 

Appendix A - Sob ols 

Aj/A^ orifice area to mainstream area ratio, 

* hr/4)/[(H/D>(S/D,j 

orifice discharge coefficient 
D orifice diameter, see Fig. 1 

Dj jet diameter ■ D 

H duct height = 10. 16 cm; see Fig. 1 

H/D duct height to orifice diameter ratio 
J momentum flux ratio, = (p j V^)/ (p^ U^) 

S spacing between centerlines of adjacent 

orifices; see Fig. 1 

S/D spacing to orifice diameter ratio 
S/H spacing to duct height ratio 
T temperature at any location in flow field 

Tj jet total temperature, ■ 300 K 

T^ mainstream total temperature, * 600 K 
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U w mainstream velocity, *15 m/sec 

Vj jet velocity 

Wj/w^ jet-to- mainstream mass flow ratio 

X distance downstieam from Jet injection loca- 

tion; see Fig. 1 

X/H dimensionless downstream distance 

Y distance above injection plane; see Fig. 1 

Y/H dimensionless vertical distance 

Z off- centerplane distance; see Fig. 1 

0 dimensionless temperature difference ratio, 

- (T^ - T)/ (T^ - Tj) 

Pj Jet density 

p x mainstream density 


establishment of a two-dimensional wall Jet flow. 
Centerline temperature difference ratio 


c, o 


EB 


*' 0 EB* 



(B3) 


where 


and 




Appendix B - Correlations 

The temperature difference ratio, 0, at any point in 
the flow field is determined from: 



The correlations for the scaling parameters in this ex- 
pression are given in this appendix. In the correlations, 
a second subscript has been appended to each of the pa- 
rameters to denote the transverse (Z) location for which 
the expression is appropriate. The scaling parameters 
are defined for planes from Z = 0 (centerplane) to 
Z * S/2 (midplane). In the conelatiom , Dj, the jet 
diameter, is used in place of D, the or.flce diameter, 
where Dj = D^cT For the range of conditions tested 
in this study, 0. 62 < < 0. 66. 

C orrelations for Predicting the Centerplane Tempera- 
tar e Profiles 


0£g is the temperature difference ratio which would re- 
sult from complete mixing of the jet and mainstream 
flows, and thus represents the asymptotic cvditlon for 
the 0 C Q decay. The decay exponent, f t varies from 
0.38 to 0.81 for the conditions tested in this study. The 
data of Ref. 6, which were for heated Jets in a confined 
crossflow' show a similar 0 C Q decay characteristic. 

In Ref. 7, a decay exponent of 0. 63 was found appro- 
priate for a single jet injected perpendicular to a semi- 
infinite crossflow. For a round Jet in a coflowing stream, 
an exponent of 0.67 would be expected. 

Minimum centerplane temperature difference ratios 

«mto.o = fl c.o( 1 - e " C T <*»> 

where 



where 


Thermal trajectory (centerline) 



•Wer 


c” =* 1. 


)/ vv 


Centerplane half-widths 


2 g ) 05 
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VL2: 


0.20 J 


.0.27. .-0.38 „ ,, 

,l t) ft) (tf » 


The exponential term is required to suppress the predic- 
ted penetration when the main flow is blocked by the 
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Note that in this investigation the jet half-widths air de- 
fined as the vertical distance from the centerline to 
where « - <« c 0 + 

Caution should be exercised in comparing the appar- 
ent profile widths directly with the half-width correlation 
relations, since the apparent widths of the profiles de- 

P endon W i/2,0- Cm. o- “ d e c.o- 

Correlations for Predicting the Off- Center-plane Varia- 
tion of file Temperature Profile Scaling Parameters 

Off-centerplane penetration 



Off-centerplane maximum temperature difference 
ratio 
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where 


d 


0.452 



Off-centerplane minimum temperature difference 
ratios and half- widths 


9 ± . 
min, z 



*1/2, z _ i/2. o 
H H 


(BIO) 


(Bll) 


Equations (BIO) and (Bll) rtpreser* "-.ajor simplifying 
assumptions which were justified by the agreement ob- 
tained between preuicted and experimental profiles. 
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Figure 5. - Variation of temperature profiles with 
distance from centerplane; X/H • 1, S/D ■ 4, 
H/D ■ 8. 
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Figure 6. - Effect of momentum flux ratio on centerplane 
temperature profiles; X/H • 1, S/D ■ 4, H/D • 12. 
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Figure 7 - El!*:t of geometric variables on centerplane 
temperature profiles at constant orifice area; 

Aj/Aoo ■ a 049. 
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Figure 8. * Effect of varying spacing on temperature 
profiles for constant orifice diamet* r ; H/D * 8. 
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Figure 9. - Effect of varying orifice diameter on centerplane 
temperature profiles at constant S/D: S/D * 2. X/H * 1. 

J - 25. 
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Figure Id - Effect of confinement by the 
opposite wall on centerplane tempera- 
ture profiles at constant S/D: S/D ■ 2. 
X/D - 4. J • 25. 
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